Lower Critical Fields of Superconducting PrFeAsOi_j, Single Crystals 
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We have studied the lower critical fields H c \ of superconducting iron oxipnictide PrFeAsOi_ H 
single crystals for H parallel and perpendicular to the afc-planes. Measurements of the local magnetic 
induction at positions straddling the sample edge by using a miniature Hall-sensor array clearly 
resolve the first flux penetration from the Meissner state. The temperature dependence of H c i 
for H || c is well scaled by the in-plane penetration depth without showing any unusual behavior, 
in contrast to previous reports. The anisotropy of penetration lengths at low temperatures is 
estimated to be ~ 2.5, which is considerably smaller than the anisotropy of the coherence lengths. 
This is indicative of multiband superconductivity in this system, in which the active band for 
superconductivity is more anisotropic. We also point out that the local induction measured at 
a position near the center of the crystal, which has been used in a number of reports for the 
determination of H c \, might seriously overestimate the obtained _ff c i-value. 

PACS numbers: 74.25.Bt, 74.25. Dw,74.25.Op,74.70.-b 



I. INTRODUCTION 

The recent discovery of high temperature superconduc- 
tivity in Fe-based compounds has attracted considerable 
interest X In this new class of compounds with a very 
low carrier density^^*^^^ superconductivity occurs 
in proximity to a magnetic instability, and unconven- 
tional pairing mechanisms mediated by magnetic fluctu- 
ations have been proposed by several groups i 10 ! 11 ! 12 One 
of the remarkable features, which is in sharp contrast to 
the high-T c cuprates, appears to be the multiband na- 
ture of superconductivity, in electron and hole pockets. 13 
Recently, a multiband effect on superconductivity has 
been reported in several compounds ! 14 ' 15 ! 16 In particular, 
the two-gap superconductivity in MgB 2 manifests itself 
in the unusual temperature- and magnetic field depen- 
dence of the anisotropy parameters in the superconduct- 
ing state i 17 ' 18 i 19 However, the crucial difference is that 
the interband coupling is very weak in MgB2 , while in Fe- 
based compounds nesting between the hole- and electron 
bands was suggested to be important for the occurrence 
of high temperature superconductivity 10 ' 11 i 12 ' 20 i 21 ' 22 i 23 
In this context, a detailed clarification of the multiband 
nature of superconductivity in the Fe-based oxypnictides 
is indispensable for the elucidation of the superconduct- 
ing properties, and especially for the pairing mechanism. 

An accurate determination of the lower critical field 
H c i is an important means to clarify not only the super- 
conducting gap symmetry, but also the multiband nature 
of superconductivity. However, the reliable measurement 
of the lower critical field is a difficult task, in particular 
when strong vortex pinning is present. We also point out 
that to date the reported values of anisotropy parameter 
strongly var y 24 i 25 i 26 i 27 spanning from 1.2 (Ref. 27) up to 



~ 20 (Ref. |24 ). which may be partly due to the effects 
of strong pinning. In this study, we use an unambiguous 
method to avoid this difficulty associated with pinning, 
by determining H c i as the field H p at which first flux 
penetration occurs from the edge of the crystal. This al- 
lows us to extract the temperature dependent values of 
the lower critical fields parallel to the c-axis (H^) and 
the afo-plane (H^), respectively as well as the anisotropy 
parameter H^/H^f in single crystals of Fe-based super- 
conductors. 



We directly determine H p by measuring the magnetic 
induction just inside and outside the edge of the single 
crystals, by using a miniature Hall-sensor array. First, 
we show that local magnetization measurements at a po- 
sition near the center of the crystal, which have been 
used by several groups for the determination of H c \, 
seriously overestimate H c \ in systems with strong pin- 
ning. Second, we find that the temperature dependence 
of H c i determined at the edge does not show any unusual 
behavior— and is well scaled by the penetration depth 
results measured on the crystal in the same batch^ 
Finally, we find that the anisotropy of the penetration 
depths 7a = A c /A Qb ~ H^/H^, where A c and A a6 
are out-of-plane and in-plane penetration depths, respec- 
tively, is much smaller than the anisotropy of the coher- 
ence lengths 75 = £ Q &/£ C = H^/H^, where £ ab and £ c 
are in- and out-of- plane coherence lengths, respectively, 
and and H^ 2 are the upper critical fields parallel and 
perpendicular to the afe-plane, respectively. This result 
provides strong evidence for the multiband nature of the 
superconductivity. 
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FIG. 1: (color online), (a) Differential magneto-optics images 
of PrFeAsOi- y #1 with field modulation SB = 0.2 mT in 
zero field. The Meissner screening occurs completely within 
narrow temperature range, (b) MO images at T = 7.1 K. 
Magnetic flux penetrates from the edge of the crystal and the 
field distribution shows the Bean critical state. 
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FIG. 2: (color online). Local magnetization loops for H\\ c, 
measured by the miniature Hall sensor located at < 10 /im 
from the edge of the crystal. 



II. EXPERIMENTAL 

Experiments have been performed on high-quality 
PrFeAsOi_ a single crystals, grown by a high-pressure 
synthesis method using a belt-type anvil apparatus 
(Riken CAP-07). Powders of PrAs, Fe, Fe 2 3 were used 
as the starting materials. PrAs was obtained by react- 
ing Pr chips and As pieces at 500°C for 10 hours, fol- 
lowed by a treatment at 850° C for 5 hours in an evac- 
uated quartz tube. The starting materials were mixed 
at nominal compositions of PrFeAsOo.6 and ground in 
an agate mortar in a glove box filled with dry nitro- 
gen gas. The mixed powders were pressed into pel- 
lets. The samples were then grown by heating the pel- 
lets in BN crucibles under a pressure of about 2 GPa 
at 1300° C for 2 hours. Platelet-like single crystals of 
dimensions up 150 x 150 x 30 /im 3 were mechanically se- 
lected from the polycrystalline pellets. The single crys- 
talline nature of the samples was checked by Laue X-ray 
diffraction.— Our crystals, whose T c (w 34 K) is lower 
than the optimum T c m 51 K of PrFeAsOi_ y ,— are in 
the underdoped regime (y ~ 0.1)^ which is close to the 
spin-density-wave order i^ 3 . The sample homogeneity was 
checked by magneto-optical (MO) imaging. MO images 
of PrFeAsOi-^ sample #1 (~ 135 x 63 x 18 ^m 3 ) are 
shown in Fig. [TJa). The crystal exhibits a nearly perfect 
Meissner state ~ 2 K below T c ; no weak links are ob- 
served, indicating a good homogeneity. At low tempera- 
tures, the magnetic field distribution is well described by 
the Bean critical state model as shown in Fig. [IJb)i2i 

The local induction near the surface of the platelet 
crystal has been measured by placing the sample on 
top of a miniature Hall-sensor array tailored in a 
GaAs/AlGaAs heterostructure^ Each Hall sensor has 
an active area of 3 x 3 /um 2 ; the center-to-center distance 



of neighboring sensors is 20 (im. The local induction at 
the edge of the crystal was detected by the miniature Hall 
sensor located at < 10 /jm from the edge. The magnetic 
field H a is applied for H\\ c and H\\ afr-plane by using 
a low-inductance 2.4 T superconducting magnet with a 
negligibly small remanent field. 

The in-plane resistivity is measured by the standard 
four-probe method under magnetic fields up to 10 T. The 
electrical contacts were attached by using the W deposi- 
tion technique in a Focused-Ion-Beam system. 



III. RESULTS AND DISCUSSION 

In Fig. [2] we show the field dependence of the "local 
magnetization", M c d gc = Mo^-^cdgc — H a , at the edge 
of the crystal, for H\\ c, measured after zero field cool- 
ing. After the initial negative slope corresponding to the 
Meissner state, vortices enter the sample and M c d gc (H a ) 
shows a large hysteresis. The shape of the magnetization 
loops (almost symmetric about the horizontal axis) in- 
dicates that the hysteresis mainly arises from bulk flux 
pinning rather than from the (Bean-Livingston) surface 
barrier i 3 ^ 

As shown in Fig. [21 the initial slope of the mag- 
netization exhibits a nearly perfect linear dependence, 
-^edge = —ctH a . Since the Hall sensor is placed on the 
top surface, with a small but non-vanishing distance be- 
tween the sensor and the crystal, the magnetic field leaks 
around the sample edge with the result that the slope 
a is slightly smaller than unity. Figure shows typical 

curves of B x l 2 = ji^ 2 (M + aHa) 1 / 2 at the edge (circles) 
and at the center (squares) of the crystal, plotted as a 
function of H a ; the external field orientation H\\ c and 
T = 22 K. The aH a -tenm is obtained by a least squares fit 
of the low- field magnetization. The first penetration field 



3 



15x10'- 



10 



-5 



H\\c 
T=22K 

• edge / ; 

■ center •* /" 

./ 

H (edge) .•/ 



crystal 

/ 



outside 



_ center edge 

W nil/ 



H (center) 
±Z l_ 



> 



4 6 
M-o^ ( mT ) 



10 



14 
12 
10 

3 8 

L= 6 

4 
2 





PrFeAsO, „#1 i 



10 



20 

r(K) 



30 



40 



FIG. 3: (color online). Typical curves of VB (left axis) at the 
edge (circles) and at the center (squares) of the crystal and 
\J Ajcdgo (right axis) plotted as a function of H a for JT|| c at 
T = 22 K, in which H a is increased after ZFC. The insets are 
schematic illustrations of the experimental setup for H || c (a) 
and if || afe-plane (b). 



FIG. 4: (color online). The temperature dependence of the 
flux penetration fields H p at the edge and the center of the 
crystal. The inset shows the temperature dependence of the 
difference between H p in the center and at the edge (left axis), 
as well as the remanent magnetization M rem (right axis). 



H p corresponds to the field ii p (edge), above which B 1 / 2 
increases almost linearly, is clearly resolved. In Fig. [31 
we show the equivalent curve, measured at the center of 
the crystal. At the center, B 1 / 2 also increases linearly, 
starting from a larger field, H p (center). 

We have measured the positional dependence of H p 
and observed that it increases with increasing distance 
from the edge. To examine whether if p (edge), i.e. H p 
measured at < 10 /im from the edge, truly corresponds to 
the field of first flux penetration at the boundary of the 
crystal, we have determined the local screening current 
density j odg c = Mo 1 (- B edge - -B outsido )/Aa; at the crystal 
boundary. Here Sedge is the local magnetic induction 
measured by the sensor just inside the edge, and -Boutside 
is the induction measured by the neighboring sensor just 
outside the edge. For fields less than the first penetra- 
tion field, jedge — f3H a is the Meissner current, which is 
simply proportional to the applied field (J3 is a constant 
determined by geometry). At H p , the screening current 
starts to deviate from linearity. Figure [3J shows the de- 
viation Aj e dge = jedge — @H a as a function of H a . As 
depicted in Fig. [3j A j c d gc again increases linearly with 
H a above if p (edge). This indicates that the if p (edge) is 
very close to the true field of first flux penetration. 

In Fig. [4l we compare the temperature dependence 
of i/p (edge) and H p (center). In the whole tempera- 
ture range, H p (center) well exceeds if p (edge). Moreover, 
H p (center) increases with decreasing T without any ten- 
dency towards saturation. In sharp contrast, H p (edge) 
saturates at low temperatures. The inset of Fig. [?] shows 
the difference between H p measured in the center and 
at the edge, AH P — H p (center) — H p (edge). AH P in- 



creases steeply with decreasing temperature. Also plot- 
ted in the inset of Fig. 0] is the remanent magnetization 
M rem (i.e. the H a = value of M edge on the decreas- 
ing field branch), measured at near the crystal center. 
This is proportional to the critical current density j c aris- 
ing from flux pinning. The temperature dependence of 
AH p is very similar to that of j c , which indicates that 
H p (center) is strongly influenced by pinning. Hence, the 
present results demonstrate that the lower critical field 
value determined by local magnetization measurements 
carried out at positions close to the crystal center, such as 
reported by several groups, is affected by vortex pinning 
effects and might be seriously overestimated] 28 ' 37 

The absolute value of H c \ is evaluated by taking into 
account the demagnetizing effect. For a platelet sample, 
H c i is given by 

H cX = H p / tanh ^0.366/0 (1) 

where a and b are the width and the thickness of the 
crystal, respectively^ In the situation where H\\ c, a — 
63 /im and b = 18 /im, while a = 18 /im and b — 63 /im 
for if 1 1 afc-plane. These values yield = 3.22H p and 
H^\ = 1.24/fp, respectively. In Fig. [51 we plot H cl as 
a function of temperature both for if || c and if || ab- 
plane. The solid line in Fig. O indicates the temperature 
dependence of the superfluid density normalized by the 
value at T = K, which is obtained from a6-plane pen- 
etration depth measurements of a sample from the same 
batch^ H^(T) is well scaled by the superfluid density, 
which is consistent with fully gapped superconductivity; 
it does not show the unusual behavior reported in Ref. |28| . 
To roughly estimate the in-plane penetration depth at 
low temperatures, we use the approximate single-band 
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FIG. 5: (color online). Lower critical fields as a function 
of temperature in PrFeAsOi-y single crystals (left axis). 
The solid line (right axis) presents the superfluid density 
Aaj(0)/A a (,(T) determined by surface impedance measure- 
ments on crystals from the same batch.— 



London formula, 



In — — 



0.5 



(2) 



where $o is the flux quantum. Using In A Q b/^ a fc + 0.5 ~ 5, 
we obtain \ a b ~ 280 nm. This value is in close corre- 
spondence with the /iSR results in slightly underdoped 
LaFeAs(0,F)^ 

Figures and (b) depict the temperature depen- 
dence of the in-plane resistivity for H [| c and H || afc- 
plane, respectively. In the inset of Fig. IHb), we dis- 
play the fields at which the resistivity is equal to 10%, 
50%, and 90% of the normal-state resistivity. For suf- 
ficiently high magnetic field, these resistance loci are 
roughly proportional to the upper critical field. In zero 
field, the resistive transition exhibits a rather sharp tran- 
sition with the transition width AT C w 2 K. By applying 
a magnetic field along the c-axis, the transition shifts 
to slightly lower temperatures and becomes broadened. 
The resistive transition curves broaden less for H\\ ab- 
plane. These results indicate that the anisotropy of the 
upper critical fields in the present system is rather large 
and that fluctuation effects play an important role for 
the transition in magnetic fields^ similar to high-T c 
cuprates4i 

Finally, Fig.[7Jshows the anisotropy of the lower critical 
fields, 7a obtained from the results in Fig. [5] Here, since 
the penetration lengths are much larger than the coher- 
ence lengths for both H\\ ab and H\\ c, the logarithmic 
term in Eq.(2) does not strongly depend on the direction 
of magnetic field. We thus assumed H^/H^ ~ A c /A a (,. 
The anisotropy 7a ~ 2.5 at very low temperature, and 
increases gradually with temperature. In Fig. the 
anisotropy of the upper critical fields 7^ is also plotted, 
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FIG. 6: (color online). Temperature dependence of the in- 
plane resistivity in PrFeAsO i- y single crystals for H || c (a) 
and if 1 1 ab- plane (b). Inset shows the temperature depen- 
dence of the upper critical fields H C 2 determined by several 
criteria that the resistivity reaches 10%, 50%, and 90% of 
the normal-state resistivity. The experimental configuration 
is also sketched. 



where 7^ is determined by the loci of 10%, 50% and 90% 
of the normal-state resistivity (see the inset of Fig.[6jb)). 
Since H C 2 increases rapidly and well exceeds 10 T just be- 
low T c for _ff || ab, plotting 7^ is restricted to a narrow 
temperature interval. In Fig. [7J we also plot the H c2 - 
anisotropy data measured on NdFeAsO0.82F0.i8 by the 
authors of Ref. [42j. These indicate that the temperature 
dependence of 7a is markedly different from that of 7^. 

According to the anisotropic Ginzburg-Landau (GL) 
equation in single-band superconductors, 7a should coin- 
cide with 7£ over the whole temperature range. There- 
fore, the large difference between these anisotropies pro- 
vides strong evidence for multiband superconductivity in 
the present system. We discuss the anisotropy parame- 
ters for the multiband superconductivity below. Accord- 
ing to GL theory, 7a and 75 at T c are given as 



7f(T c ) = 7A(T c ) 



(3) 



where (• • •) denotes the average over the Fermi surface, v a 
and v c are the Fermi velocities parallel and perpendicu- 
lar to the a6-plane, respectively. 43 ' 44 f2 represents the gap 
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FIG. 7: (color online). Normalized temperature depen- 
dence of the anisotropies of H c \ (-jx, closed circles) and H C 2 
(75, closed squares) in PrFeAsO i_ a single crystals. The 
anisotropy of H C 2 in NdFeAsO0.82F0.i8 (7^, open squares) 
measured by Y. Jia et al^ is also plotted. The dashed line 
is a guide to the eye. 



isotropic scattering^ which corresponds to 7a ~ 4. This 
value is close to the observed value. The fact that 7^ 
well exceeds 7a indicates that the active band for super- 
conductivity is more anisotropic than the passive band. 
According to band structure calculations, there are five 
relevant bands in LaFeAsOi^Fa;. Among them, one of 
the three hole bands near the T point and the electron 
bands near the M point are two-dimensional and cylin- 
drical. The other two hole bands near the T point have 
more dispersion along the c axis,— although the shape of 
these Fermi surfaces is sensitive to the position of the As 
atom with respect to the Fe plane, which in turn depends 
on the rare earths Our results implying that the active 
band is more anisotropic is in good correspondence with 
the view that the nesting between the cylindrical hole 
and electron Fermi surfaces is essential for superconduc- 
tivity. This is expected to make these two-dimensional 
bands the active ones, with a large gap, and the other 
more three-dimensional bands passive ones with smaller 
gaps. 



IV. SUMMARY 



anisotropy ((^ 2 ) = 1), which is related to the pair poten- 
tial V(v, v') = V fl(v)Q(v'). At T = K, the anisotropy 
of the penetration depths is 

75(0 = ]f . (4) 

The gap anisotropy does not enter 7a (0), while 7^ at 
T = K is mainly determined by the gap anisotropy of 
the active band responsible for superconductivity. Thus 
the gradual reduction of 7a with decreasing temperature 
can be accounted for by considering that the contribu- 
tion of the gap anisotropy diminished at low tempera- 
tures. This also implies that the superfluid density along 
the c-axis X 2 (0)/X 2 (T) has steeper temperature depen- 
dence than that in the plane X 2 b (0)/X 2 b (T). A pro- 
nounced discrepancy between 7^ and 7a provides strong 
evidence for the multiband nature of superconductivity in 
PrFeAsOi-y, with different gap values in different bands. 
We note that similar differences between 75 (T) and 
7 a(T), as well as \ 2 C (0) / \ 2 C (T) and \ 2 ab (0) / \ 2 ab (T) , have 
been reported in the two-gap superconductor MgB 2 i 18 ' 19 
We also note that angle-resolved photoemission spec- 
troscopy (ARPES), 45 Andreev reflection,— and pene- 
tration depth 4 ^ measurements on (Ki_ :E Ba 2 ,)Fe2As2 and 
and penetration deptbJS, studies of LnFeAs(0,F) 
(Ln = Pr, Sm) have suggested multiband superconduc- 
tivity with two gap values in Fe-based oxypnictides. 

Band structure calculations for LaFeAsOi-^Fz yield 
an anisotropy of the resistivity of approximately 15 for 



In summary, we have measured the lower critical field 
H c i in PrFeAsOi_j, single crystals for H\\ c and Jf || ab- 
plane by utilizing an array of miniature Hall sensor. Con- 
ventional methods using a single micro-Hall probe placed 
on the center of the crystal might overestimate H c \ due 
to strong flux pinning. H c \ measured by the sensor lo- 
cated very near to the edge of the crystal shows satu- 
rating behavior at low temperatures, which is consistent 
with the previous reports on the penetration depth mea- 
surements. The anisotropy of H c \ slightly decreases with 
decreasing temperature and is indicative of multiband su- 
perconductivity in PrFeAsOi_j,, in which the active band 
for superconductivity is more anisotropic. 
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